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We investigate the hydrodynamic low-frequency oscillations of highly degenerate Fermi gases 
trapped in anisotropic harmonic potentials. Despite the lack of an obvious spatial symmetry the 
wave-equation turns out to be separable in elliptical coordinates, similar to a corresponding result 
K-j , established earlier for Bose-condensates. This result is used to give the analytical solution of the 

■ anisotropic wave equation for the hydrodynamic modes. 

<+H 1 I. INTRODUCTION 

o 

After the spectacular realization of Bosc-Einstcin condensation in evaporatively cooled magnetically trapped alkali 
vapors half a decade ago ||] the experimental realization of a degenerate (i.e. mean particle distance smaller than 
the de Broglie wavelength) Fermi gases like 40 if Q and 6 Li || has also been achieved. One of the most useful 
experimental tools to probe the properties of the Bose-condensates has been the excitation of their collective modes 
in the collisionless regime . It has turned out that within the dynamic Thomas- Fermi approximation these modes 
could be calculated analytically by solving the Gross-Pitaevskii equation, linearized around the static Thomas-Fermi 
£j ■ solution. This was done first |(^,|7j for isotropic parabolic traps, where the symmetries are sufficient for integrability, 
l — 1 \ then for anisotropic axially symmetric parabolic traps |^,|[@] , where an additional hidden symmetry could be identified 
[gj , and finally also for completely anisotropic parabolic traps |lO[ ] , where the resulting wave-equation could be related 
to various integrable dynamical systems and the Lame equation. 
. In view of these successful experimental and theoretical developments for the Bose-systems, it seems to be of 
considerable interest to apply similar methods also to degenerate Fermi systems On the one hand these 

systems seem to be very different from their bosonic counterparts. Indeed, barring the possibility of Cooper pairing 
at extremely low temperatures, condensates may not appear there. However, a degenerate Fermi gas may appear, 
instead. Therefore a great similarity between both kinds of systems remains, because in the strongly degenerate case 
the particles move in the trap-potential and strong mean-fields coming from their effective interactions, either with the 
condensate in the bosonic case, or via the Pauli exclusion-principle with the degenerate gas, in the case of fermions. 

In the present paper our goal is to analyze the wave-equation of the collective modes of the degenerate Fermi gas 
in a completely anisotropic (i.e. with three different trap-frequencies) parabolic trap. More symmetric special cases 
^j" 1 , of this problem have already been analyzed in the literature. 

For the special case of an isotropic parabolic trap the wave-equation for hydrodynamic modes in trapped Fermi 
gases was solved by Bruun and Clark |l2[ both in the high- and low-temperature limit. In the high-temperature limit 
the effects of quantum statistics become irrelevant and one obtains the same results as for Bose gases, which were 
considered earlier, for isotropic traps, and even for some modes in axially symmetric traps by Griffin et al. jl4|. Using 
a sum-rule approach Vichi and Stringari Jll| analyzed interaction effects on the collective oscillations of Fermi gases 
trapped in parabolic potentials of rotational or axial symmetry. Amoruso et al. fl3f| also investigated the collective 
excitations of trapped degenerate Fermi-gases both in the hydrodynamic and the collisionless regime for isotropic and, 
for some dipolar and quadrupolar modes, also for axially symmetric parabolic traps. New in the present work is the 
treatment of the completely anisotropic case, where we follow the lines of our earlier treatment of the bosonic case 
[fL0|. The paper is organized as follows. 
In section |0| the hydrodynamic equations are formulated and the wave-equation for small density perturbations in 



the limit T —> is derived. In section III the wave-equation is solved in the axially symmetric case by separation in 



cylindrical-elliptic coordinates. We also construct the conserved operator commuting with the wave-operator, whose 



presence is responsible for the integrability. In section IV we proceed to the triaxially anisotropic case and separate 



the wave-operator in three-dimensional elliptic coordinates, exhibiting two conserved operators commuting with the 
wave-operator. Analytical expressions for the mode-frequencies and -eigenfunctions are given. Finally, we summarize 
our conclusions. 
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II. EQUATIONS OF MOTION 



The hydro-dynamic equations of motion of an ideal gas in local thermodynamic equilibrium with energy-density 
e(x, t) = |-P(x, t), where P denotes the pressure, are given by ]l5| , [T^ | 

du VP 1 - , . 

p— = p-VV(x) 2.1 

at mm 

— = --(v-uP + pu-VV{x)) (2.2) 

with d/dt = d/dt + u-V. 

Here p is the total density of the gas of fermions, which we assume to be present in (2F + 1) different hyperfmc- 
states, which all see the same trapping potential and between which s-wave collisions are not suppressed by the 
Pauli-principle. We shall not consider here additional hydrodynamic modes connected with the internal degrees of 
freedom of the fermions, like spin- waves. 

The trapping potential is assumed to have the form 

^)^(^+^+^). (2.4) 

The number-density po (x) in equilibrium is determined by the Fcrmi-Dirac distribution in the semi-classical approxi- 
mation 

fo(x,p) = , 1 , (2.5) 

e P(^+ v ( s )-^ + 1 

where f3 = l/(fcsT) and p is the chemical potential. fo(x,p) reduces in the low-temperature limit to 

f (x,p) = e(^-^--V(x)), (2.6) 
\ 2m / 

where Q(x) is the Heaviside jump-function. Thus, summing over the hyperfine-states a we get 



Pol 

which , for T — > 0, reduces to 



(*)=£/7^/o(^1> (^) 



(2F+ 1) 
6tt 2 



Z7T7 

-t(p-v(x)) 

n 



3/2 



Q(p-V(x)), (2.8) 



where &(■) is the Heaviside jump-function. The degenerate Fermi-gas has a sharp boundary at p = V(x) which forms 
a Fermi-ellipsoid for parabolic trapping potentials in complete analogy to the Thomas-Fermi ellipsoids formed by large 
condensates in such traps. Integrating over the whole cloud of particles we get 

with the geometric mean trap-frequency 

uj = (u x u y u z ) 1/3 . (2.10) 
The pressure Po(x) in equilibrium satisfies the hydrostatic condition 

VPo(x) = - Po {x)VV{x), (2.11) 
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from which 



P (x) = -( t M-V(x))p (x) 
5 



follows. The same result follows, for T — > 0, from 



2 2(2^+1) r d 3 fc ft 2 fc 2 
P (z) = r o(x) = y —ffaM). 



We can eliminate /x — V{x) from eqs.( 2.12] ) and ( 2.1§| ) to obtain the equation of state 



Po 



ih 2 



6tt 2 



5 to V 2F + 1 



2/3 



5/3 



Linearizing eqs.(2.1), ( |2.2| ), and (2.3) around equilibrium, with p = po + Sp, P = Pq + (5Pq we obtain 



= --V5P - Sp-VV(x) 
at to m 



(2.12) 



(2-13) 



(2.14) 



(2.15) 



<95P _ 5 - 
~df ~ ~3 



P (z)u 



-p (x)u-VV(x) 



(2.16) 



= -V 



p (x)u 



(2.17) 



Eliminating i5P and dp from these equations using the hydrostatic equilibrium condition relating Pq(x) and V(x) and 
the property po(x) = p (V(x)) we obtain fll4| , |l2"[ 



= f (v • ff ) - f (w) v • a- If (*• w) 

c/t 3to £>n V / 3m V / to V / 



9i 2 3to po ^ I 3m 
Using in addition eq.( 2.12| ) we obtain the wave-equation 



d 2 u _ 1^ 
di 2 m 



_( M -F(x))V-u - (V-V(x) 



(2.18) 



(2.19) 



Excluding time-independent flows, the gradient-form of the right-hand side permits the introduction of a velocity- 
potential ip 



u = —\7ip 

in terms of which the wave-equation becomes 

dtp=-£-Qi- V{x)) VV - ~ (W(f)) ■ V*>. 
3m to V / 



d 2 ip 
dt 2 



From the continuity equation we can derive, with u = —\7ip and dSp/dt =/= 0, 

3to Po(x) dip 
P ~ 2 p-V(x)~dt' 

Similarly, from the energy-balance equation it follows that 

SP = mpo(x)^-, 



which implies the relation 5P = |(p — V(x))Sp, and shows, together with eqs.( |2.12 ) and ( p. 18 ), that 



sp = dP " {po) s P . 

opo 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



(2.24) 
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In the case of an isotropic trap the potential V(x) becomes V(x) = ima^r 2 . Rescaling lengths by the Thomas-Fermi 



radius ttf — y and separating the time-dependence by <p(x, t) — e lu)t ip(x) the wave-equation (2.21) in spherical 
coordinates r, <p, 6 simplifies to 



uj 2 ip(r, (f>, 9) = u 2 



l (l-r 2 ) fd 2 d 

3 r 2 \dr dr e " 



dr 



(p(r, 0, 0). 



Here 77V 2 , i is the angular part of the Laplacian. The separation ansatz ip — r l F n j(r 2 )Yi ri 
harmonics Y\ m , combined with the ansatz of a power-series of finite order for F n ; 



(2.25) 

with the spherical 



F nl (r 2 )=^C k r 2k , 



k=0 



yields for the term of highest power 2k = 2n in 



,2(+2n 



1 d 9 d 1(1 + 1) d 

r 2 + r — 

3 dr dr 3 dr 



J+2n 



from which the mode-frequencies 



2 4 , / 2 , „ 31 

to = -Un n +nl + 2n H 

3 u V 4 



(2.26) 



(2.27) 



(2.28) 



follow, as first derived by Bruun and Clark (jT2| and Amoruso et al. 0|. The polynomials F n i(x) are given by the 
hypergeometric functions |12| 



F n i (x) — const ■ 2-F1 (— n,l + 2 + n;l + -;x). 



(2.29) 



In the following sections we consider the solution of eq.( |2.21 ) for axially symmetric and for triaxially anisotropic 
parabolic traps. 



III. SOLUTION FOR AXIALLY SYMMETRIC TRAPS 



In order to solve eq.(2.21) for axially symmetric trapping potentials 

V(x) = \^mu\ (x 2 + y 2 ) + ^mu) 2 z 2 
it is useful to introduce cylindrical-elliptic coordinates 

\/ x 2 + y 2 = (7 \f (1 + £ 2 )(1 — V 2 )! z — a i 1 l^ 4> = arctan—, 

x 

The focal length a of the Thomas-Fermi ellipsoid is 



I2u / 1 1 



(3.1) 



(3.2) 



(3.3) 



its excentricity is e = yT — Lo\jio 2 . The expressions as given apply for the case u>± < ui z . The case w_L > u>\\ can be 

treated similarly with results which can also be obtained by the analytical continuation of the results given here, see 
H. The wave-equation can be transformed to these coordinates. Separating the time-dependence via 



it takes the form 

(l_ e 2„ e 2^2 )(1 _ e 2 +£ 2 ?/ 2 ) J 2 

(l-e 2 )(e + V 2 ) 



<p(2, t) = e- iwt (p(x) 



r) 2 fl 2 f) f) 



(3.4) 



drj 2 <9£ dr] \ 1 — rf 1 + £ 
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2e 2 (l + £ 2 ) d_ 2e 2 (l-7? 2 ) 8_ 



1 



l-e 2 -e 2 £ 2S <9£ l-e 2 + eV'a77 
In order to solve this equation we make the ansatz 

which permits the separation of variables and leads to the two equations 



1 - e 2 - e 2 ^ 2 1 - e 2 + e 2 r) 2 



2w 2 I 

7T^ = ° 



2 I",. ,9n d 2 d w 2 



2e 2 (l+£ 2 ) t d 



l-e 2 -e 2 e 2 rfC l-e 2 -e 2 ^ 2 



(3.5) 



(3.6) 



(3.7) 



.,, d 2 d 
2i] 



drj 2 drj 1 — rf 



2e 2 (l-r/ 2 ) d 2w 2 /w 2 



1 - e 2 + e 2 r] 2 drj 1 - e 2 + e 2 ?7 2 



^(r?) = -B<p v (n). 



(3.8) 



Here -B is the separation constant. It is the eigenvalue of a conserved operator i3, commuting with the wave-operator 
G, 



[B,G]=0. 



(3.9) 



It can be constructed from eqs.(3.7) and (3Js) by eliminating the eigenvalue ui 2 of G from both equations and replacing 
the eigenvalue m of the angular-momentum L around the z-axis by the operator L z = —id/d(f). We obtain in elliptic 
coordinates 



2/ 9 , ,J , 9. 9 „ 9 , 9 2/1 1 

-3 (1 - £ + eV) ^ (1 " ^ " 2£ (1 - " } % + 3 (iT^ - W 



(3.10) 



Converted to the original coordinates B takes the simple form 



B= 2 - 
3 



(f.v) (. 



x- V + 4 - 



2/i / a 



a 2 



2^ a 2 



mw^ \ dx 2 dy 2 ) mtu 2 dz 2 



(3.11) 



The operator B can be freed from any dependence on the trap geometry ||l7| (apart from the parabolic form of the 
trap-potential, which is essential for its existence) by rescaling the coordinates x — > x 



x = xj A 



2//. 



V = y/\ 



2/i 



2 ' 



z = zL 



2(i 



(3.12) 



whence it becomes 



B 



2 r 



(*■*) 



x- V + 4 - V 



(3.13) 



which is rotationally invariant in the rescaled coordinates. Apart from a trivial factor 2/3 the operator B differs from 
its bosonic counterpart found in || merely by a term proportional to x • V. 

The eigenvalue-problem of B can be easily solved simultaneously with that of L z by the ansatz, in cylindrical 
coordinates p, cj>, z, 



= e irn<f> p \m\ zaF (m, a ) ^2 ^ 



(3.14) 



where a = 0, 1 determines the parity in z and F^™ 1 ^ is a polynomial of finite order n 



5 



n k 



'p J Z 



(3.15) 



fc=0 j=0 

The condition that n remains finite determines B as 



B = - (2n + a + |m|) (2n + a + \m\ + 4) . 



(3.16) 



Since B and L z both commute with the wave-operator G they also have simultaneous eigenfunctions. Moreover, as 
we shall now show, the polynomial eigenfunctions of B with non-vanishing eigenvalue B ^ satisfy the requirement 
of particle-number conservation 



d 3 xSp(x, t) = 0. 



By eqs.(2.22), pj) and (RJ), eq.(3.17) for w ^ is equivalent to 



• / -•• i i x2 y 2 z2X 1 ' 

iuj I a X ( 1 tt x r 

fl z fl z c z 



y>(x) = 



or, in rescaled space-coordinates and for uj ^ 0, 

d 3 x(l-x 2 ) 1/2 lf (Z) = 0. 



(3.17) 



(3.18) 



(3.19) 



Mu ltiply ing with B ^ and inserting the operator to replace its eigenvalue we can evaluate the left-hand side of 
eq.( 3.19| ) by partial integration 



B / d 3 x(l - x 2 y/ 2 <p= / d 3 x(l-x 2 ) 1/2 



V • xx ■ V - V + x ■ V 



2 3 x (v(i 



= / d 3 x 



( /' ; .;;(V(l-f 2 ) 
; i _ l 



•2x1/2 



ii • V - V ) ip + I d 3 x (1 - x A ) H " x ■ V(f 



X 



-X ' 



X 



■Vtp + d 3 x(l - i 2 ) 1/2 x- Vt, 



= 0. 



(3.20) 



We conclude that only the eigenfunctions with a — = \m\ — n must be excluded, because they correspond to B = 
and uj = 0. 

We must also guarantee momentum conservation of the modes, which means in a harmonic oscillator potential 
that the total momentum P(t) = J d 3 xpo(x)u(x,t) must satisfy P{ + uifPi — 0, where we denoted the three main 
frequencies u/ x , oj v , lo z by u>i, i = 1,2,3 and the Cartesian components of P by Pj. Using the wave-equation ( 2. IS ) 
and our results for po(x), Pq(x) we can calculate 



Of 2 



d 3 xp u = +- d 3 xVi(P V ■ u) 



dh 



V-(p uViV)-Vi(p o u-fV 



(3.21) 



which can give only surface-terms at the boundary of the degenerate Fermi-gas. These vanish for polynomial modes 
u = —Vip because Pq ~ (p — V(x)) 5 ^ 2 and po ~ (p — V(x)) 1 ^ 2 vanish at the boundary p — V(x). Hence 



0L 

dt 2 



d xpou = 0. 



(3.22) 



The wave-equation (2.21) leads to a recursion relation for the coefficients C^.™'"' 1 in a straight-forward manner. The 



coefficients corresponding to the highest total power k — n satisfy a tridiagonal homogeneous set of equations 



__ o 2 

j>=0 ^ 



(3.23) 



G 



with the matrix 



M jf = 5 jtj , 



4j(j + |m|)+6j + 3H | 4(rc-j + f)(rc-j + f + 1) 



? {n- J + - + l 



a 1 

n — i H 1 — 

J 2 2 



(3.24) 



+^-'-i--o(i + 1 )0' + 1 + H) 



where we denote the radial and axial half-axis of the Fermi-ellipsoid as a = 2[i j rauj\ and c = y/2/x/ma/| , respectively. 
For any given pair of integers n, \m\ and parity a = 0, 1 this finite homogeneous set of (n + 1) equations is easily 
solved, in general numerically, and determines (n + 1) eigenfrequencies ^™'' m '' Q ' labeled by the additional integer 
quantum number j with < j < n. 

Simple solutions are obtained for n = 0. We have already shown that a = m = must be excluded if n — 0, due 
to particle-number conservation. Therefore, the smallest frequencies are 



(0,0,1) 

= W|| 



,(0,1,0) _ 



which are the frequencies of the Kohn- modes |16) . More generally 



(0,|m|,0) 



with the (unnormalized) number-density modes 



5p{x) = consul I- 4 - ^e im V m| z a - 



1/2 



The eigenfrequencies for n = 1 are also easily obtained and given by 

W G) H,0 ~ Vf ( 5 ^ (|m| + 1} + 4w t ± x/t 2 !^! + 5 )M + H^M Z (H + 4) ^ } ) 

upj mW = ^1 (5wi(|m| + 1) + 9wf ± ^(2|m| + 5) 2 ^ + - ^ )j V ' . 

The number-density modes take the form 

The relation between the coefficients C 1 ™' a ' > and c[™' a ^ for the n — 1 modes can be obtained from eq.( 3.23 ) 



C 



, v , , 3w 2 - (7|m| + 10)w 2 - wjfa 



10 



2^1(1 + 2a) 



(3.25) 



(3.26) 



(3.27) 



(3.28) 



(3.29) 



(3.30) 



and Cq™'"" 1 is then most easily obtained from the recursion-relation of the eigen-value equation for B, which yields 



rM _ 3 a 2 (|m| + l)C { ;r ] + c 2 (l/2 + a)Cftr ,a) 
00 " 2 3+|m| + a 

The results for the eigenfrequencies ( |3.2g| ) for the special case \m\ = 0, a = have already been obtained in Jl 



(3.31) 
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IV. TRIAXIALLY ANISOTROPIC CASE 



Let us now turn to the case of a parabolic potential without any further symmetries. We shall analyze this case using 
the procedure applied to Bose-condensates in our previous paper (To). Ellipsoidal coordinates X,p, v are introduced 
as functions of x, y, z as the three roots of |0| with respect to p 



V 



c 2 + p 



1 



er + p tr + p 

where a> b> c are the semi-axis' of the Fermi-ellipsoid. We order the roots according to 

- a 2 < v < -b 2 < fi < -c 2 < A < 0. 
In order to avoid confusion the chemical potential is denoted as p from now on. It is useful to define 

F(p) = (a 2 +p) (b 2 + p) (c 2 +p) 
The wave-equation can then be written as 

( F(X) 1 



m 9,9 o 9 , 4A/i^ 

-aWwV = - 

p (A — p)(p — v)(y — A) 



(fi - v) 



2 d 2 

1™W ' V A 



-F'(X) 



8_ 

8X 



+ cyclic. 



It can be separated like in the bosonic case |10(] with the result 



^F(A)J^+(F(A) + iAF'(A) 



3_ 

8X 



if : 



ttllo 2 a 2 b 2 c 2 AX BX 2 
— z — + — + — — 

p 4 4 4 



(4.1) 

(4.2) 
(4.3) 

(4.4) 
(4.5) 



and two cyclic equations in p, v. A and B are separation constants. They are the eigenvalues of operators A and B 
which commute with the wave-operator G. The existence of these conserved operators is the reason for the solvability 
of the wave-equation des pite the absence of any apparent symmetries of the anisotropic potential. In order to construct 
A and B we replace in ( |4.5| ) A and B by their operators and solve these three equations for Aip and Bip. Afterwards 
we transform back to Cartesian coordinates. We obtain in this way after some calculation 



A = I { [(b 2 + c 2 )(x 2 - a 2 ) + a\y 2 + z 2 )] ^ + 2a\jz ^ ' ' ^ 9 



,- 4(6 + c )x— + cyclic 
oyoz ox 



(4.6) 



which differs from the bosonic case only by the trivial prefactor 2/3 and the factor 4 instead of 3 in the last term. For 
the operator B we obtain 



B 



$ ■ V ( f • V + 4^ 



d 



2 1 



a 



B 2 B 2 

By 2 8z 2 



8x 2 



(4.7) 



which is a simple generalization of the result ( 3.11 ), and reduces again to ( 3.15 ) after rescaling x,y,z by a, 6, c, 
respectively. The eigenfunctions of B and its eigenvalues are therefore still given by polynomials of finite total order 
Tit in x, y, z and 



B = -n t (n t + 4), 



n t = 0,1,2, 



(4.8) 



The eigenfunctions of the wave-operator can now be obtained as in |10| using the method for the solution of the Lame 
equation (see 18|). The polynomial ansatz (with n t = In + a + j3 + 7; a = 0, 1; j3 = 0, 1; 7 = 0, 1) 



ip = x a y l3 z' r Y]_ ( Pi(x) 



(4.9) 



with undetermined coefficients 6s. where 



tpi[x) 



b 2 



c z + 1 



(4.10) 
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inserted in the wave-equation Gip — ~uj 2 ip, yields after some rearrangement 



— U) if 



2a 2/3 2 7 

H 



n , 



x 2 y 2 z 2 \ ^ dip(x) 



b 2 



5 



4 4(2a+l) 4(2/3+1) 4(2 7 +l) 16 
^ + 3 a 2 + e t + 3 6 2 + 6>i + 3 c 2 + 0] + T 



n 



(4.11) 



Here we define 



dip(x) 
dipi(x) 



— T "J»7 



n ^( f )- 



3 = 1 
W#0 



The ansatz for ip is a solution if the second term on the left-hand side of eq.( 4.11 ) vanishes. This is the case if, for 
i = 1, 2, ...rt 



6 4a + 2 4/3 + 2 4 7 + 2 
& + a 2 +6» 4 + 6 2 + (9i + c 2 + 9i 



E 



= 0. 



J=1 '» "J 



It can be shown that all the roots of ( 4.12 ) are real [|10|. The eigenfrequencies are then given by 



2 2/1 / a P 7 



n 

2 E^ 



(4.12) 



(4.13) 



We consider the cases n = and n = 1. For n = we find the 8 frequencies, for a, 0, 7 = 0, 1 



2fi f a (3 7 



with 



(4.14) 



(4.15) 



The case a = /3 = 7 = with uj = must be excluded however. The same frequencies and modes are obtained i n the 
bosonic case. For n > 1 the mode-frequencies in the bosonic and fermionic case differ, however. For n — 1 eq.( 4.12| ) 
reduces to a cubic equation for #i = 9 



(4(a + /? + 7) + 10) (6 1 + a 2 )(0 + 6 2 )(0 + c 2 ) - [(4a + 2) a 2 (0 + b 2 ) (9 + c 2 ) + cyclic] = 
and the mode-frequencies for the eight possible choices of the a = /3 = 7 are 



2Jx ( a (3 7 2 



(4.16) 



(4.17) 



In the general case n > 1 the 6^, i = 1,2, ...n can be determined as the equilibrium of the n-particle potential 



2 + i)Eh« 2 



l + \)±in\b 2 + t \-(l + \)±in\c 2 + e t 



(4.18) 



n n 



E S> 

2 — 1 J— 2 + 1 



<) 



which differs only in the prefactor of the first term on the right-hand side from the bosonic result (cf. eq. (4.24) 
of jLO)). To find the minima of the potential ( 4.18 ) is a problem in the 2-dimensional electrostatics (because of the 
logarithmic potential) of n+4 point-charges on a line: The equilibrium positions 9i of the n-charges of unit strength 
between the 4 fixed charges of strength 3/4 a t = 0, (f + |) at = -a 2 , (f + \) at 6 = -b 2 , and + i) at Q = -c 2 
determine the eigenfrequencies via eq.( 4.13| ). The solutions can be classified by the numbers ni,n2,n3 of charges 
between — c 2 and 0, — b 2 and c 2 , and —a 2 and — b 2 , respectively, which are the three quantum numbers labeling the 
eigenvalues and eigenfunctions. We shall not present numerical results here, since these depend on the ratios a/6, a/c 
as an input. For a further discussion we can refer to En], however, where the bosonic case is analyzed further by 
numerically solving the equivalent electrostatic problem for a specific experimentally realized trap-geometry. 



V. CONCLUSIONS 



In the present paper we have analyzed the hydrodynamic equations of an ideal and highly degenerate Fermi-gas 
trapped in a generally anisotropic harmonic potential. Despite the absence of any obvious spatial symmetries in 
the triaxially anisotropic case, or of insufficient spatial symmetry in the axially symmetric case, the wave-equation 
turns out to be separable just like in the bosonic case. In fact, from the point of view of dynamical systems the 
characteristics of the wave-equation (i.e. the rays described by the Hamilton-Jacobi equation obtained in the cikonal 
limit [^0|) are the same for the bosonic and the fermionic case; i.e. the wave-equations for the two cases can be 
considered as different 'quantizations' of the same 'classical' dynamical system, which was analyzed in [JlO|, and which 
turned out to be integrable. 

The frequencies of the lowest-lying modes the n = modes (4.14), (4.15) (but not the mode-functions for Sp(x)) 
turn out to be the same in the bosonic and the fermionic case. The reason for this coincidence is the fact that the 
velocity-potential (p of these modes satisfies V 2 </? = 0, which implies the incompressibility condition V • u = . Fo r 
incompressional modes the only terms depending on the quantum-statistics drop out from the wave-equation ( 2.18| ) 
for the velocity-field, as was already noted by Griffin et al. |Q. However, the frequencies of the higher-lying modes 
differ, in general, between the Bose- and the Fermi-systems. In the axially symmetric case we have given all the 
mode-frequencies with at most 1 node in radial or axial direction, and with arbitrary axial angular momentum. In 
the completely anisotropic case useful explicit results for compressional modes with V • u ^ can only be given 
in numerical form. However, as a central result of the present work we have reduced the numerical task to a very 
straight-forward one, namely to the 2-dimensional electrostatics of n movable point-charges of unit strength on a line 
in the field of 4 fixed point charges of given strength on the same line. This establishes a strong analogy between the 
fermionic case and the bosonic case, which has been reduced to the same problem (with slightly different strengths of 
the fixed charges) in [|l(| . As discussed there this mapping of the eigenvalue-problem permits a complete classification 
of the eigenvalues and the eigenfunctions by the three integer quantum- numbers, which give the number of movable 
charges between each pair of nearest neighbors of the 4 fixed charges. 

The physical realization requires the simultaneous trapping of several hyperfine states, in order to allow for the 
interaction of the different fermionic species by collisions, which are necessary to maintain local thermodynamic equi- 
librium. Because the collision-rates in degenerate Fermi-gases are suppressed by a Fermi-blocking factor proportional 
to (T/Tp) 2 compared to the classical collision rates [[l9| it is necessary to use atomic species with particularly large 
positive or negative s-wave scattering-lengths or to produce the latter by tuning the system close to a Feshbach 
resonance. 

In view of the recent rapid experimental progress in this field Q , || it may be hoped that a measurement of these 
collective modes may be possible in the near future. 
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